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Abstract 

Background and Aim: Defective apoptosis is a hallmark of cancer development and progression. Death receptors (DR4, 
FAS) and their ligands (TRAIL, FASL) are thought to mediate the major extrinsic apoptotic pathway in the cell. SNPs in these 
genes may lead to defective apoptosis. Hence, the present study aimed to investigate the association of functional SNPs of 
DR4 (rs20575, rs20576 and rs6557634), FAS (rs2234767) and FASL (rs7631 10) with gallbladder cancer (GBC) risk. 

Methods: This case-control study included 400 GBC and 246 healthy controls (HC). Genotyping was carried out by Taqman 
genotyping assays. Statistical analysis was performed by using SPSS ver16. Meta-analysis was performed using 
Comprehensive Meta-analysis software (Version 2.0, BIOSTAT, Englewood, NJ) to systematically summarize the possible 
association of SNP with cancer risk. Functional prediction of these variants was carried out using Bioinformatics tools (FAST- 
SNP, F-SNP). False discovery rate (FDR test) was used in multiple comparisons. 

Results: The DR4 C rs20575 A rs20576 A r5655763 4, G rs20575 A rs20576 G rs6557634 and G rs20575 C rs20576 G rs655753 4 haplotypes conferred two- 
fold increased risk for GBC. Among these, the DR4 C rs20 5 75 A rs20 5 76 A rs6 55 76 34 haplotype emerged as main factor influencing 
GBC susceptibility as the risk was not modulated by gender or gallstone stratification. Our meta-analysis results showed 
significant association of DR4 rs6557634 with overall cancer risk, Gl cancers as well as in Caucasians. We didn't find any 
association of FAS and FASL SNPs with GBC susceptibility. 

Conclusions: The DR4 haplotype C rs20 5 7 5A rs2 o5 7 6A rs6 55 76 34 represents an important factor accounting the patients 
susceptibility to GBC probably due to decreased apoptosis. However, additional well-designed studies with larger sample 
size focusing on different ethnicities are required to further validate the results. 
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Introduction 

Gall bladder carcinoma (GBC) is the most common biliary tree 
cancer and fifth most common gastrointestinal malignancy, with 
striking geographical, age, race, and gender-related differences in 
its incidence [1]. Incidence of GBC in north India is 7.4 per 10' 
for females and 3.6 per 10 for males in Delhi [2]. The higher 
incidence of GBC in females has been attributed in part to 
hormonal factors, including estrogen and progesterone receptors. 
GBC is characterized by rapid progression and a very high 
mortality rate [3,4]. Diagnosis of over 90% of GBC patients at an 
inoperable stage with serious invasion and metastasis to other 
organs is ascribed to absence of specific signs or symptoms [5] . 
Inadequate information and poor understanding regarding 
molecular basis of GBC [6] necessitate the establishment of 
important inroads in understanding the pathophysiology, early 
diagnosis and inter-individuals variability in susceptibility of GBC. 



TRAIL-R 1 (TNF-related apoptosis-inducing ligand receptor 1)/ 
DR4 and Fas/Apo-l/CD95 are members of tumor necrosis factor 
receptors which link exogenous stimuli via transmembrane surface 
receptors to the intracellular signaling machinery that mediates 
and executes the death signal. These pathways are one of the 
major extrinsic apoptotic signaling pathway [7] and defects in 
death receptor signaling can confer resistance to apoptosis. The 
role of apoptosis in tumorigenesis has been well-documented [8] 
and resistance to apoptosis is believed to be a hallmark of cancer 
[9] . Thus TRAIL/TRAIL-R and Fas/FasL system seems to play a 
crucial role in the pathogenesis of cancer. DR4 is thought to be an 
attractive candidate tumor suppressor gene [10]. DR4 and FAS 
mutations resulting loss or reduction of apoptotic function have 
been described in different human cancers [11,12]. Aberrant 
expression of FAS and/ or FASL has been detected in many 
human cancers and appears to be a feature of the malignant 
phenotype [13,14]. 
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Since single-nucleotide polymorphisms (SNPs) are the most 
common forms of human genetic variation and may affect cancer 
risk by influencing individual susceptibility [15], SNPs in DR4, 
FAS and FASL gene that impair the apoptotic signals represent 
the plausible and promising etiologic pathways modifying the 
GBG penetrance and merit to be tested as risk factors of 
gallbladder cancer susceptibility. The rs20575, located in exon 4, 
immediately 3 ' to the main receptor ligand interface regions of the 
DR4 protein, results in the substitution of an arginine for 
threonine. It was suggested to influence the receptor ligand 
binding resulting deficient apoptotic signaling [11]. The rs20576 
replaces glutamate by alanine which leads to structural change 
within region of TRAIL/DR4 complex formation. This may result 
in insufficient interaction with TRAIL and obstructed induction of 
caspase-8 dependent apoptosis, causing longer survival rate of 
tumor cells [16]. The FAS -1377 G to A substitution has been 
shown to decrease the expression of FAS protein [17] while the C 
allele of FASLG -844 C>T has been shown to increase the 
expression of the FASL protein [18]. Recently, many genetic 
studies have investigated the role of DR4, FAS and FASL 
polymorphisms in the etiology of various cancers [16,19-30] etc. 
Nevertheless, the results of these studies remain inconclusive 
because some cancer types were in positively associated with these 
polymorphisms, while others were not. Therefore, we have carried 
out the present case-control study with the goal to find whether 
polymorphism of these apoptotic proteins modulates GBC 
susceptibility in North Indian population. 

Materials and Methods 

Ethics statement 

The study protocol was sanctioned by the ethical committee of 
Sanjay Gandhi Post Graduate Institute of Medical Sciences 
(SGPGIMS) Lucknow (India) and King George Medical Univer- 
sity (KGMU), Lucknow. The study subjects were registered 
according to the norms of the World's Association Declaration of 
Helsinki. Demographic characteristics data, such as gender, age, 
occupation history, tobacco history and other lifestyle factors, as 
well as clinical and pathological investigations were collected in a 
custom-designed questionnaire by conducting a personal interview 
with all study participants. Staging of cancer was done according 
to the AJCC-TNM classification system [31]. 

Subjects 

A total of 400 histo-pathologically confirmed gallbladder cancer 
patients from north Indian population were included in the 
present study. All the patients were enrolled from the Department 
of Gastroenterology of Sanjay Gandhi Postgraduate Institute of 
Medical Sciences, a tertiary care hospital in Lucknow, Uttar 
Pradesh, India and Department of Surgical Oncology, KGMU, 
Lucknow, Uttar Pradesh, India belonging to same ethnicity. For 
controls, 246 healthy individuals (age and sex matched with similar 
ethnicity to cases) unrelated to patients and to each other were 
recruited in the study from volunteers visiting the hospital for a 
routine checkup or health awareness camps and hospital 
employees from same region to ensure similar ethnicity. Samples 
were collected from June 2006 to Dec 2012. 

The family history and general questionnaire of all the patients 
recruited for the study were taken to assure their ethnicity. The 
healthy controls had no evidence of any personal history of 
gallstone, gallbladder cancer or other malignant conditions or 
chronic disease. The inclusion and exclusion criteria for the 
recruitment of the Patients and controls were as previously 
described [32]. The participation rate was 100%, and blood 



samples were available for all subjects. 3 ml of blood was collected 
in EDTA vial and stored at — 70°C. Informed and written 
consents were taken from all participants recruited in the study 
when interviewing for the demographic details and blood sample 
collection. 

DNA isolation 

Genomic DNA was extracted from peripheral blood using the 
standard salting out method [33]. The quality and quantity of 
DNA were checked by using Nanodrop spectrophotometer 
(Thermo Fisher Scientific/Nanodrop Products, Wilmington, 
Delaware, USA). 

SNP selection 

We have selected three coding and functional polymorphisms of 
DR4 gene (C>G/Thr209Arg/rs20575, A>C/Glu228Ala/ 
rs20576 and G >A/Argl41His/rs6557634 which are the most 
commonly studied polymorphism. In addition, promoter poly- 
morphism of FAS (-1377G/A, rs2234767) and FASL (- 844T/C, 
rs763110) believed to be involved in regulation of transcription 
were also included in the study. 

Genotyping of DR4, FAS and FASL polymorphisms 

Taqman allelic discrimination assay on ABI 7500 Real Time 
PCR system using 96-well plates were used for genotyping of 
SNPs: rs20575, rs20576, rs2234767 and rs7631 10. Pre-designed 
assay primers and probes were purchased from Applied Biosys- 
tems (Foster City, CA, USA). A negative control (wells containing 
no DNA) was included in all plates. Amplification was detected by 
using ABI PRISM 7500 sequence detection software. SNP 
rs6557634 was analyzed using AS-PCR (Allele-specific Polymerase 
Chain Reaction). Details of the primers and cycle conditions for 
this SNP were taken from earlier published work [28] . Positive and 
negative controls were used in genotyping assay, and 5% of the 
samples were randomly selected and run in duplicates with 100% 
concordance. The results were consistent with no discrepancy in 
genotyping. 

Statistical analysis 

Sample size was estimated via Quanta version 1.1.1 [34] using 
HapMap (http:/ /hapmap. ncbi.nlm.nih.gov/) minor allele fre- 
quency data to attain 80% power. Hardy Weinberg Equilibrium 
in the control population was tested by employing Goodness of fit 
X 2 test. 

All statistical analysis was performed using SPSS statistical 
analysis software, version 17.0 (SPSS, Chicago, IL, USA). 
Haplotypes estimation and linkage disequilibrium analysis was 
conducted using the SNPstat Software [35]. FAST-SNP (http:// 
fastsnp.ibms.sinica.edu.tw) and F-SNP (http://compbio.cs. 
queensu.ca/F-SNP) was used for in-silico study. 

Descriptive statistics of patients and controls were presented as 
the mean and standard deviations (SDs) for continuous measures, 
while frequencies and percentages were used for categorical 
measures. The differences in demographic variables and genotype 
distributions of the polymorphisms between cases and controls 
were compared by using Chi- square analysis or two-sided Fisher's 
exact test. The most common homozygous genotype was taken as 
a reference to calculate the risk estimates for codominant, 
dominant, and recessive genetic models. Binary logistic regression 
analysis was used for all analysis variables to estimate odds ratio 
(OR) and 95% confidence interval (CI) adjusted for age and 
gender to estimate the risk of gallbladder cancer with the 
polymorphisms. A two-tailed p-value of less than 0.05 was 
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considered a statistically significant result. To keep meaningful 
results and maintain a low rate of false-positive findings, FDR was 
applied in case of multiple statistical testing [36]. 

Meta-analysis 

Literature was searched out in PubMed, Scopus and Embase, 
using the keywords: T>R4/rs20575/rs20576/rs6557634', 'poly- 
morphism', 'cancer'. The search was updated extending all the 
published papers until July 2013. Reference lists of key studies and 
reviews were also hand-screened for additional potentially eligible 
studies. All the studies were included in the present meta-analysis if 
they met the following criteria: a) original papers, b). human- 
association studies, c). exploring the association between any of the 
three selected SNPs of DR4 gene and cancer susceptibility, d). 
case-control studies, e). papers with crude Odd's ratio and 95% 
confidence interval or sufficient data to calculate overall OR at 
95% CI). Unpublished findings, reports without control popula- 
tion or genotype frequency data, duplication of previous publica- 
tions and non-English articles were excluded. Information on the 
several data such as: first author's name, publication date, country, 
cancer type, ethnicity, number of patients/control included in the 
study, genotyping methods etc. was collected for each study. 
Ethnic groups were categorized as Asian, European and Mixed. 

In meta-analysis, the strength of association between DR4 SNPs 
and risk of cancer was assessed by calculating the crude odds ratio 
(ORs) and corresponding 95% confidence interval (CI). Analyses 
were weighted by trial size. The pooled ORs were estimated for 
allele contrast (Variant vs. wild allele), heterozygote versus wild 
genotype, the homozygote variant genotype versus wild homozy- 
gote and also in dominant model. Statistical heterogeneity was 
evaluated using the 0_ statistic among the studies and it was 
considered as significant when p<0.05 [37]. Only random effects 
model (Der Simonian Laird) was used since it is more conservative. 
Publication bias was investigated with Begg's funnel plot. Funnel 
plot asymmetry was further assessed by the method of Egger's 
linear regression test with p<0.05 being considered statistically 
significant [38]. Moreover, sensitivity analysis was also performed, 
excluding studies whose allele frequencies in controls exhibited 
significant deviation from the Hardy-Weinberg equilibrium 
(HWE), given that the deviation may denote bias [39]. The 
significance of overall odds ratio (OR) was determined by the Z- 
test. All the p-values were two sided and all the statistical tests were 
performed using the Comprehensive Meta-analysis software 
(Version 2.0, BIOSTAT, Englewood, NJ). 

Results 

Population Characteristics 

Table 1 presented the demographic characteristics of GBC 
patients with respect to their age and gender matched controls. 
The mean age of 400 GBC and 246 controls are 52.65± 10.45 
year and 47.75± 10.65 year, respectively. There was no statistically 
significant difference among the mean age of controls and cases. 
About 95% of the GBC patients were in advanced stages of cancer 
(stage III and stage IV) and 31% of the GBC patients were tobacco 
users (either in some form smoking, chewing, or both). About 
50.0% of GBC patients were associated with gallstones. All cancer 
patients were incident cases and none of the controls had a family 
history of cancer. Among GBC, 37% of the cases had early age of 
onset, i.e. <50 years. 

Allelic distribution of studied polymorphisms in controls 

The distributions of DR4, FAS and FASL polymorphisms 
(rs20575, rs20576, rs6557634, rs2234767 and rs763110) are 



Table 1. Demographic profile of the study subjects. 





GBC 


HC 


Variables 


N (%) 


N (%) 


Subjects 


400 


246 


Gender 


Male 


122 (30.5) 


83 (33.7) 


Female 


278 (69.5) 


163 (66.3) 


Age ±SD (year) 


52.65 ±10.45 


47.75 ±10.65 


Gallstone 


No 


200 (50.0 ) 


246 (100) 


Yes 


200 (50.0) 




Tobacco 


No 


273 (68.9) 




Yes 


123 (31.1) 




Stages 


0, 1 


No (0) 




II 21 (5.25) 


III 199(49.75) 


IV 


1 80 (45.0) 





HC - Healthy Control, GBC - Gallbladder cancer. 
doi:1 0.1 371 /journal.pone.0090264.t001 



shown in Table 2. The observed genotype frequencies of the 
two studied polymorphisms in HC were in accordance with 
Hardy-Weinberg equilibrium (p>0.05). 

Association of DR4, FAS and FASL polymorphisms with 
GBC 

The risk of GBC in relation to each of the SNPs studied in DR4, 
FAS and FASL genes are shown in Table 2. On comparing the 
genotype frequency distribution of DR4 SNPs (rs20575, 
rs20576and rs763110) in GBC patients with that of controls, we 
found a significant association of rs20576 with GBC risk at 
heterogenotype (AC, pcorr = 0.043, OR= 1.74), variant allele (C, 
pcorr = 0.033, OR = 1 .7 1) as well as in dominant model (AC+CC, 
pcorr = 0.035, OR =1.78). Similarly, heterozygous carriers of 
DR4 SNP rs6557634 also showed a significant association with 
increased gallbladder cancer risk (GA, pcorr = 0.035, OR= 1.61). 
This association also persisted at allele level (A, pcorr = 0.030, 
OR = 1 .45) as well as in dominant model (GA+AA, pcorr = 0.028, 
OR= 1.66). In contrast, there were no significant differences in 
genotype frequencies between gallbladder cancer cases and 
controls for DR4 rs20575. 

In haplotype analysis, the DR4 haplotype comprising the major 
alleles was taken as reference and the difference in the frequencies 
of haplotypes between patients and controls were tested using chi- 
square test (Table 3). The haplotype analyses of the three 
studied SNPs of DR4 revealed that the frequencies of 

Crs20575A r s20576A r5 65 5 7 6 34; G rs 20575A rs2 0576G r s6557634 an d 

G rs 20575C rs 2o.'i76Gr S 6557634 haplotypes were significantly higher in 
GBC as compared to controls (pcorr = <0. 0005, OR = 2.76; 
pcorr = 0.0244, OR = 2.09 and pcorr = 0.0335, OR = 2.80, 
respectively). These haplotypes were found to confer two-fold 
increased risk for GBC. The global haplotype analysis also 
demonstrated statistically significant differences (pcorr = 0.0015) 
between GBC cases and controls based on the distribution pattern 
of the all haplotypes (Table 3). 
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Table 2. Overall frequency distribution of DR4, FAS and FASL polymorphism in GBC and HC. 





Genotype/Allele 


HC N (%) 


GBC N {%) 


OR (95%CI) 


p-value 


FDR pcorr 


DR4 rs20575 genotypes/alle 


es (age and gender 


adjusted) 








GG 


128 (52.0) 


183 (45.8) 


Reference 






CG 


101 (41.1) 


181 (45.2) 


1.05 (0.69-1.59) 


0.827 




CC 


1 7 (6.9) 


36 (9.0) 


1.94 (0.90-4.15) 


0.090 




P-trend 








0.104 




CG+CC 


118 (48.0) 


217 (54.2) 


1.16 (0.78-1.72) 


0.472 




GG+CG vs. CC 


229 (93.1) 


364 (91.0) 


1.90 (0.91-3.97) 


0.090 




G 


357 (72.6) 


547 (68.4) 


Reference 






C 


135 (27.4) 


253 (31.6) 


1.22 (0.90-1.67) 


0.198 




DR4 rs20576 genotypes/alle 


es (age and gender 


adjusted) 








AA 


210 (85.4) 


300 (75.0) 


Reference 






AC 


33 (13.4) 


86 (21.5) 


1.74 (1.03-2.95) 


0.039 


0.043 


CC 


3 (1.2) 


14 (3.5) 


2.19 (0.47-10.19) 


0.320 




P-trend 








0.001 




AC+CC 


36 (14.6) 


100 (25.0) 


1.78 (1.07-2.95) 


0.026 


0.035 


AA+AC vs. CC 


243 (98.8) 


386 (96.5) 


1.97 (0.42-9.17) 


0.390 




A 


453 (92.1) 


686 (85.8) 


Reference 






C 


39 (7.9) 


114 (14.2) 


1.71 (1.08-2.71) 


0.023 


0.033 


DR4 rs6557634 genotypes/a 


lleles (age and gene 


er adjusted) 








GG 


126 (51.2) 


1 74 (43.5) 


Reference 






GA 


107 (43.5) 


1 92 (48.0) 


1.61 (1.06-2.44) 


0.025 


0.035 


AA 


13 (5.3) 


34 (8.5) 


2.05 (0.90-4.70) 


0.088 




P-trend 








0.029 




GA+AA 


120 (48.8) 


226 (56.5) 


1.66 (1.11-2.48) 


0.013 


0.028 


GG+GA vs. AA 


233 (94.7) 


366 (91.5) 


1.63 (0.73-3.62) 


0.234 




G 


359 (73.0) 


540 (67.5) 


Reference 






A 


133 (27.0) 


260 (32.5) 


1.45 (1.07-1.98) 


0.018 


0.030 


FASL rs7631 10 genotypes/a 


leles (age and gender adjusted) 








TT 


90 (36.6) 


145 (36.2) 


Reference 






CT 


119 (48.4) 


180 (45.0) 


0.97 (0.63-1.50) 


0.898 




CC 


37 (15.0) 


75 (18.8) 


1.14 (0.64-2.03) 


0.667 




P-trend 








0.048 




CT+TT 


1 56 (63.4) 


255 (63.8) 


1.01 (0.67-1.53) 


0.947 




TT+CT vs. CC 


209 (85.0) 


325 (81.2) 


1.15 (0.68-1.96) 


0.598 




T 


299 (60.8) 


470 (58.8) 


Reference 






C 


193 (39.2) 


330 (41.2) 


1.05 (0.79-1.40) 


0.739 




FAS rs2234767 genotypes/a 


leles (age and gend 


?r adjusted) 








GG 


153 (62.2) 


245 (61.2) 


Reference 


- 


- 


AG 


86 (35.0) 


1 36 (34.0) 


0.95 (0.62-1 .44) 


0.802 




AA 


7 (2.8) 


1 9 (4.8) 


1 .68 (0.56-5.07) 


0.356 




P-trend 








0.537 




AG+AA 


93 (37.8) 


1 55 (38.8) 


1.00 (0.67-1.50) 


0.997 




GG+AG vs. AA 


239 (97.2) 


381 (95.2) 


1.72 (0.56-5.11) 


0.333 




G 


392 (79.7) 


626 (78.2) 


Reference 






A 


100 (20.3) 


174 (21.8) 


1.06 (0.75-1.49) 


0.753 





HC - Healthy Control, GBC - Gallbladder cancer, OR - Odds Ratio, CI - Confidence Interval. 
Significant Values are given in bold. 
doi:1 0.1 371 /journal.pone.0090264.t002 
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Table 3. Frequency distribution of DR4 haplotypes in GBC patients and HC. 



Hctplotypes 


Frequency 






OR (95%CI) 




HC 246 (%) 


GBC 400 (%} 






Crs20575A rS 20576G r s6557634 


0.5917 


0.4316 




Reference 


Crs20575A r s20576A r s6557634 


0.119 


0.2095 


<0.0001 a 


2.76 (1.71 - 4.47) 


G r s20575Ars20576G r s6557634 


0.0947 


0.1427 


0.0084 b 


2.09 (1.21 - 3.62) 


G r s20575A r s20576Ars6557634 


0.1154 


0.0737 


0.62 


0.86 (0.47 - 1.57) 


G r s20575^rs20576G rS 6557634 


0.0306 


0.0719 


0.022 c 


2.80 (1.16-6.76) 


G r s20575^rs20576A r s6557634 


0.0337 


0.028 


0.23 


1.75 (0.71-4.32) 


Crs2057sCrs20576Grs6557634 


0.0127 


0.0288 


0.31 


1.83 (0.57-5.92) 


Global haplotype association p-value: 0.00038 d 



HC - Healthy control, GBC - Gallbladder cancer, OR - Odds Ratio, CI - Confidence Interval. 
Significant Values are given in bold, FDR pcorr = 0.0005", 0.0244 b , 0.0335°, 0.0015 d . 
doi:1 0.1 371 /joumal.pone.0090264.t003 



However, no significant differences were observed for, FAS 
rs2234767 and FASLrs763110 SNP in any groups both at 
genotypic and allelic levels as well as for dominant and recessive 
models (Table 2). 

Association of DR4, FAS and FASL polymorphisms with 
GBC after Gender stratification 

Stratification of subjects on the basis of gender showed that the 

C rs 20575Ars20576Ars6557634 an d G rs 2O575^rs20576G r s6557634 haplo- 

types of DR4 were significandy associated with a two-fold 
increased risk of GBC in females (pcorr = 0, OR = 2.60 and 
pcorr = 0.043, OR=2.10, respectively]. However, the C rs 20575 
A rs 20576Ars6557634 haplotype was found to confer increased risk of 
GBC in males (pcorr = 0.0277, OR = 3.55) as compared to 
controls (Table 4). The global haplotype association was also 
significant only in female GBC patients (pcorr = 0.0245) as 
compared to controls. 

No association was observed between FASL SNP rs7631 10 and 
GBC susceptibility in both males and females. However, the 
variant allele of FAS SNP rs2234767 was found to confer 
significantly protective effect for GBC in males (A; pcorr = 0.043, 
OR = 0.51, Table SI in File SI). 

Association of DR4, FAS and FASL polymorphisms with 
GBC after gallstone stratification 

Since GBC is often associated with gallstones, we stratified 
patients on the basis of presence/ absence of gallstone to see the 
combined effect of gallstone and DR4, FAS and FASL variants on 
GBC risk. A significant difference was found for DR4 C rs 20575 

A rs 20576A rs 6557634j an d G rs 20575C rs 20576A rs 6557634 haplotypes ill 

GBC with stones as compared to controls (pcorr = <0. 0005, 
OR = 4.04 and pcorr = 0.0277, OR = 3.78, respectively). Similar- 

tyi C rs 20575Ars20576Ars6557634> G rs 20575A. s 20576G r s6557634 and 

G rs 20575C rs 20576G rS 6557634 haplotypes of DR4 gene were found to 
significantly higher in GBC without gallstone as compared with 
control (pcorr = 0.0282, OR=1.97; pcorr = 0.0154, OR=2.35 
and pcorr = 0.0434, OR=2.81, respectively). The global haplo- 
type association in GBC patients with and without stone was also 
significant (pcorr = <0.0005 and 0.0277, respectively; Table 5). 

Table S2 in File SI shows the distribution of SNPs after 
stratifying the subjects on the basis of gallstone. The frequencies of 
FAS and FASL genotypes did not differ significantly between 



control and GBC with or without stones, showing lack of 
association (Table S2 in File SI). 

Meta-analysis for DR4 rs20575 

There were fourteen studies investigating the association of 
DR4rs20575 with cancer risk [10,19,21,22,25,28,40-46] but we 
have excluded two studies of Kazuya Kuraoka 2005 [43] and 
Veronica Fernandez 2004 [44] because these studies either 
deviated from HWE (p<0.05) [43] or used tissue samples [44] 
for polymorphism analysis. We combined the results of our study 
with the rest 12 eligible studies [10,19,21,22,25,28,40-42,45,46] to 
calculate pooled OR for 4731 cases/5029 control and found no 
significant overall association with risk of cancer [C vs. G: 
OR= 1.05 (0.97-1.12), p = 0.222; CC vs. GG: OR= 1.07 (0.94- 
1.22), p = 0.291; CG vs. GG: OR= 1.08 (0.92-1.27), 0.352; 
CC+CG vs. GG: OR= 1.09 (0.94-1.26), p = 0.268; Table 6 and 
Figure SI in File SI]. On subgroup analysis, no significant 
association was found amongst any groups after ethnicity or cancer 
site based subgroupings. 

Meta-analysis for DR4 rs20576 

For rs20576, 8 studies were found eligible [10,16,19,22,28,40- 
42]. Wolf et al, 2006 [16] conducted the association in CLL, 
MCL, bladder cancer, prostate cancer and HNSCC. These studies 
were taken as individual study because of different type of cancer; 
however we have excluded the studies involving HNSCC and 
prostate cancer as tissue samples were used in analysis. Thus, a 
total of 1 1 studies (Table S3 in File SI) including our study, having 
a total of 2589 cases/ 3222 controls (400 patients/246 controls 
from present study) were included. The meta-analysis results 
indicated heterogeneity, hence random model was used for meta- 
analysis [phet = <0.05, Figure S2 in File SI]. OveraU, this SNP 
did not show any association with cancer risk [C vs. A: OR= 1.23 
(0.95-1.56), p = 0.112; CC vs. AA: OR = 1.48 (0.85-2.57), 
p = 0.162; AC vs. AA: OR= 1.16 (0.87-1.53), p = 0.310; CC+AC 
vs. AA= 1.21 (0.91-1.62), p = 0.192]. In subgroup analysis, also, 
rs20576 did not show any significant association with cancer risk. 

Meta-analysis for DR4 rs6557634 

For rs6557634, only three studies were found [19,28,42] and on 
combining the results of these studies with our study (total 832 
patients/723 control), overall a marginal significant association 
was found for cancer risk with A vs. G [OR= 1.22 (1.05-1.42), 
p = 0.011], AA vs. GG [OR= 1.52 (1.06-2.19), p = 0.023] and in 



PLOS ONE | www.plosone.org 



5 



February 2014 | Volume 9 | Issue 2 | e90264 



Death Receptor Variants and Gallbladder Cancer 



Table 4. Haplotype Analysis of DR4 gene in GBC and HC after gender stratification. 



Haplotypes 


Haplotype Analysis of DR4 in 


male 


Haplotype Analysis of DR4 in females 


Frequency 






Frequency 








HC (83) 


GBC (122)p- value 


OR (95%CI) 


HC (163) 


GBC (278) p- value 


i— \ n men/ f ' * 
OR (95%CI) 


Crs20575A r s20576G r s6557634 


0.5826 


0.4181 




Reference 


0.6008 


0.4366 




Reference 


Crs20575Ars20576A r s6557634 


0.1309 


0.2351 


0.01 3 a 


3.55 (1.32-9.54) 


0.1079 


0.1984 


0.0000 b 


2.60 (1.49 -4.52) 


G r s20575A r s20576Grs6557634 


0.1122 


0.1435 


0.069 


2.44 (0.94-6.34) 


0.084 


0.1435 


0.035 c 


2.10 (1.06 - 4.17) 


G r s20575Ars20576A r s6557634 


0.108 


0.0681 


0.74 


0.81 (0.22-2.90) 


0.1215 


0.0758 


0.066 


0.86 (0.43-1.70) 


G r s20575^rs20576G rS 6557634 


0.016 


0.0447 


0.092 


9.10 (0.71-117.02) 


0.0392 


0.0803 


0.073 


2.31 (0.93-5.78) 


G r s20575^rs20576A r s6557634 


0.0349 


0.0265 


0.34 


0.57(0.07-4.89) 


0.0314 


0.0313 


0.12 


2.39 (0.80-7.19) 


Crs2057sCrs20576G r s6557634 


0 


0.0495 


0.61 


5.57 (0.17-186.44) 


0.0153 


0.0231 


0.59 


1.43 (0.38-5.36) 


Crs2057sCrs20576A r s6557634 


0.0154 


0.0145 


0.36 


4.78 (0.17-136.99) 


NA 


0.011 


0.0000 





Global haplotype association p-value= 0.11 



Global haplotype association p-value =0.0092 



GBC - Gallbladder cancer, HC - Healthy controls, OR - Odds Ratio, CI - Confidence Interval 
Significant Values are given in bold, FDR pcorr = 0.0277 a , 0 b , 0.0430 c , 0.0245 d . 
doi:1 0.1 371 /journal.pone.0090264.t004 



dominant model [AA+AG vs. GG: OR= 1.29 (1.05-1.60), 
p = 0.017, Figure S3 in File SI]. In subgroup analysis these 
association were found limited to only Caucasians and for GI 
cancer (Table 6). 

Bias diagnostics 

Funnel plot shape for the pooled analysis did not reveal any 
asymmetry in all comparison models. Egger's test also did not 
suggest any evidence of publication bias for rs20575 [CC vs. GG: 
t= 1.28, p = 0.228; CG vs. GG: t = 2.12, p = 0.057; CG+CC vs. 
GG: t= 1.89, p = 0.086; and C vs. G allele: t = 1.13, p = 0.281, 
Figure S4 in File SI], for rs20576 [CC vs. AA: t= 1.49, p = 0.170; 
AC vs. AA: t = 0.66, p = 0.525; AC+CC vs. AA: t = 0.90, 
p = 0.392; C vs. A allele: t = 1.08, p = 0.308, Figure S5 in File 
SI] and for rs6557634 [AA vs. GG: t = 0.91, p = 0.460; AG vs. 
GG: t = 0.47, p = 0.687; AG+AA vs. GG: t = 1.87, p = 0.202; A vs. 
G aUele: t = 2.82; p = 0.106, Figure S6 in File SI]. 



In-silico analysis 

Table S4 in File SI show our in-silico analysis for the five 
studied SNPs. The rs20575 was found to be involved in splicing 
regulation, though the effect was tolerated and benign (FS 
score = 0.330). Similarly, rs6557634 and was found to be involved 
in splicing, possibly damaging and benign (FS score = 0.284). The 
rs20576 was also found to be possible damaging but tolerated. The 
rs763110 was involved in transcriptional regulation while 
rs2234767 was found benign. 

Discussion 

TRAIL-Rl and Fas death receptors form a subgroup of the 
tumor necrosis factor receptors (TNF-R) superfamily having a 
conserved cytoplasmic signaling module (—80 amino acid) termed 
the death domain DD which is essential for apoptosis induction 
[47]. Binding of these receptors to their cognate ligands (i.e., 
TRAIL and FASL, respectively) lead to receptor aggregation and 
recruitment of adaptor proteins, and caspases, forming the Death- 



Table 5. Haplotype analysis of DR4 gene in GBC (with stone and without stone) and HC. 



Haplotype Analysis of DR4 in GBC without stone and 

HC Haplotype Analysis of DR4 in GBC with stone and HC 

Frequency Frequency 



Haplotypes HC (246) GBC (200) p- value OR (95%CI) HC (246) GBC (200) p- value OR (95%CI) 



^rs2057sArs20576G r s6557634 


0.5917 


0.4271 




Reference 


0.5917 


0.4399 




Reference 


Crs2057sA r s20576Ars6557634 


0.119 


0.1978 


0.015 a 


1.97 (1.14 -3.41) 


0.119 


0.2165 


<0.0001 e 


4.04 (2.23 -7.32) 


G r s20575A r s20576Grs6557634 


0.0947 


0.1666 


0.0048 b 


2.35 (1.30-4.23) 


0.0947 


0.115 


0.24 


1.58 (0.74- 3.35) 


G r s20575Ars20576A r s6557634 


0.1154 


0.0659 


0.32 


0.69 (0.34 - 1 .42) 


0.1154 


0.086 


0.86 


1.07 (0.50 - 2.30) 


G r s20575^rs20576G rS 6557634 


0.0306 


0.0737 


0.042 c 


2.81 (1.04-7.61) 


0.0306 


0.0695 


0.075 


2.61 (0.91-7.46) 


G r s20575^rs20576A r s6557634 


0.0337 


0.0187 


0.94 


0.95 (0.27-3.33) 


0.0337 


0.0369 


0.01 3 f 


3.78 (1.32-10.79) 


Crs2057sCrs20576G r s6557634 


0.0127 


0.0375 


0.21 


2.26 (0.64-7.99) 


0.0127 


0.0205 


0.52 


1 .60 (0.38-6.67) 




Global haplotype association p-value: 0.016 d 


Global haplotype association p-value: <0.0001 9 



GBC-Gallbladder cancer, HC-Healthy controls, OR -Odds Ratio, Cl-Confidence Interval 

Significant Values are given in bold, FDR pcorr = 0.0282 a , 0.01 54 b , 0.0434 c , 0.0284 d , 0.0005 e , 0.0277 f , 0.0005 s . 
doi:1 0.1 371 /journal.pone.0090264.t005 
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Table 6. Meta-analysis Results. 



DR4 rs20575 



Variables 


N 


Case/ 
Control 


C vs G allele 






CC vs GG 






CG vs GG 






CC+CG vs GG 












OR (95% CI) p 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


All 


13 


4731/ 
5029 


1 .05 (0.97- 
1.12) 


0.222 


0.203/ 
23.764 


1 .07 (0.94- 
1.22) 


0.291 


0.363/ 
8.328 


1.08 (0.92-1.27) 


0.352 


0.830/ 
0.000 


1.09 (0.94-1.26) 


0.268 


0.011/ 
53.862 


Ethnicity 


Caucassian 


9 


4304/ 
4578 


1 .06 (0.98- 
1.16) 


0.153 


0.095/ 
40.807 


1.11 (0.95- 
1.30) 


0.192 


0.200/ 
27.484 


1.07 (0.89-1.29) 


0.446 


0.002/ 
67.245 


1.10 (0.93-1.30) 


0.282 


0.004/ 
64.711 


Mixed 


3 


292/282 


1 .04 (0.82- 
1.31) 


0.764 


0.594/ 
0.000 


1 .08 (0.67- 
1.72) 


0.764 


0.586/ 
0.000 


1.31 (0.82-2.10) 


0.258 


0.295/ 
18.173 


1.22 (0.82-1.81) 


0.323 


0.418/0.000 


Asian 


1 


135/169 


0.87 (0.63- 
1.20) 


0.380 


1 .000/ 
0.000 


0.79 (0.43- 
1.44) 


0.435 


1 .000/ 
0.000 


0.84 (0.50-1.41) 


0.511 


1.000/ 
0.000 


0.82 (0.51-1.32) 


0.413 


1.000/0.000 


Cancer type 


Gl Cancer 


6 


1862/ 
1821 


1.12 (0.98- 
1.28) 


0.086 


0.116/ 
43.318 


1 .26 (0.98- 
1.62) 


0.071 


0.237/ 
26.340 


1.15 (0.85-1.56) 


0.354 


0.002/ 
73.820 


1.18 (0.90-1.54) 


0.235 


0.004/ 
70.920 


Other Cancer 


7 


2869/ 
3208 


0.99 (0.92- 
1.06) 


0.728 


0.796/ 
0.000 


0.98 (0.85- 
1.14) 


0.786 


0.804/ 
0.000 


0.98 (0.84-1.14) 


0.761 


0.250/ 
23.434 


0.96 (0.86-1 .08) 


0.516 


0.436/0.000 


DR4 rs20576 






























Variables 


N 


Case/ 
Control 


C vs A allele 






CC vs AA 






AC vs AA 






CC+AC vs AA 












OR (95% CI) 


P 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


All 


11 


2589/ 
3222 


1.23 (0.95- 
1.59) 


0.112 


0.000/ 
81.444 


1 .48 (0.85- 
2.57) 


0.162 


0.014/ 
54.807 


1.16 (0.87-1.53) 


0.310 


0.000/ 
76.133 


1.21 (0.91-1.62) 


0.192 


0.000/ 
79.793 


Ethnicity 


Caucassian 


8 


2297/ 
2940 


1.32 (0.97- 
1.80) 


0.076 


0.000/ 
86.343 


1 .60 (0.87- 
2.95) 


0.131 


0.007/ 
63.824 


1.19 (0.85-1.67) 


0.300 


0.000/ 
81 .886 


1.28 (0.91-1.82) 


0.162 


0.000/ 
85.004 


Mixed 


3 


292/282 


0.98 (0.69- 
1.41) 


0.930 


0.342/ 
6.577 


0.98 (0.20- 
4.75) 


0.977 


0.255/ 
26.719 


1.02 (0.61-1.71) 


0.931 


0.213/ 
35.309 


1.00 (0.64-1.55) 


0.994 


0.289/ 
1 9.400 


Cancer type 


Gl Cancer 


5 


1653/ 
1574 


1 .24 (0.84- 
1.84) 


0.285 


0.000/ 
87.721 


1 .43 (0.70- 
2.91) 


0.327 


0.030/ 
62.643 


1.18 (0.79-1.76) 


0.430 


0.000/ 
81 .897 


1.23 (0.79-1.92) 


0.358 


0.000/ 
86.120 


Other Cancer 


6 


936/1648 


1.23 (0.82- 
1.86) 


0.319 


0.001/ 
75.590 


1.91 (0.57- 
6.41) 


0.295 


0.049/ 
54.975 


1.14 (0.70-1.84) 


0.608 


0.002/ 
73.815 


1.20 (0.76-1.91) 


0.436 


0.001/ 
74.698 


DR4 rs6557634 


Variables 


N 


Case/ 
Control 


A vs G allele 






AA vs GG 






AG vs GG 






AA+AG vs GG 












OR (95% CI) 


P 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


OR (95% CI) 


P 


Ph/I 2 


All 


4 


832/723 


1.22(1.05- 
1.42) 


0.011 


0.604/ 
0.000 


1.52(1.06-2.19)0.023 


0.390/ 
0.367 


1.24 (0.99-1.54) 


0.059 


0.830/ 
0.000 


1.29 (1.05-1.60) 


0.017 


0.909/ 
0.000 


Ethnicity 


Caucassian 


2 


600/471 


1.31 (1.08- 
1.58) 


0.005 


0.923/ 
0.000 


2.03 (1.23- 
3.34) 


0.006 


0.771/ 
0.000 


1.26 (0.97-1.62) 


0.079 


0.748/ 
0.000 


1.34 (1.05-1.71) 


0.019 


0.875/ 
0.000 


Mixed 


2 


232/252 


1.07 (0.83- 
1.37) 


0.623 


0.658/ 
0.000 


1.11 (0.66- 
1.88) 


0.692 


0.588/ 
0.000 


1.18 (0.76-1.83) 


0.464 


0.397/ 
0.000 


1.16 (0.77-1.77) 


0.480 


0.668/ 
0.000 


Cancer type 


Gl Cancer 


2 


600/471 


1.31 (1.08- 
1.58) 


0.005 


0.923/ 
0.000 


2.03 (1.23- 
3.34) 


0.006 


0.771/ 
0.000 


1.26 (0.97-1.62) 


0.079 


0.748/ 
0.000 


1.34 (1.05-1.71) 


0.019 


0.875/ 
0.000 


Other Cancer 


2 


232/252 


1 .07 (0.83- 
1.37) 


0.623 


0.658/ 
0.000 


1.11 (0.66- 
1.88) 


0.692 


0.588/ 
0.000 


1.18 (0.76-1.83) 


0.464 


0.397/ 
0.000 


1.16 (0.77-1.77) 


0.480 


0.668/ 
0.000 



Significant associations are shown in bold, Ph- p-value of Q test for heterogeneity, OR -Odds Ratio, Cl-Confidence Interval. 
doi:1 0.1 371 /journal.pone.0090264.t006 



Inducing Signaling Complex (DISC). This, in turn, activate and 
releases caspase 8 and 10 for the triggering of apoptosis [48]. 
Normal variations within the sequence of apoptotic genes are 
suggested to lead suboptimal apoptotic capacity, finally increasing 
cancer risk [49]. Despite many investigations, the associations 



between DR4, FAS, and FASL polymorphisms and the risk of 
human cancers remain inconsistent between different studies. 
Therefore, in this study, we tested the association of coding 
polymorphisms of the death receptor 4 (DR4), FAS, and its legend 
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(FASL) with GBC susceptibility. In addition a meta-analysis was 
also performed to find out the overall cancer risk. 

The rs20575, was suggested to result deficient apoptotic 
signaling [1 1]. However, we did not find any association of DR4 
rs20575 polymorphism with the GBC risk. Our results are in 
agreement with other studies showing null association of rs20575 
with lung cancer [28], breast cancer [10,25,42,45], bladder cancer 
[19], and gastric cancer [43,50]. In our meta-analysis also, we 
didn't find any association with cancer risk. A previous meta- 
analysis by Chen et al., 2009 also indicated only marginal 
association of rs20575 with overall cancer susceptibility [51]. Our 
in-silico polyphen analysis has also showed rs20575 to be benign. 

For DR4 rs20576, we have observed an increased risk of GBC 
at hetero-genotype level, allelic level and in dominant model. Our 
in-silico polyphen study also showed rs20576 to be possibly 
damaging. The rs20576 variant was found to exhibit an enhanced 
risk to have CLL, MCL, HNSCC and bladder cancer as well as an 
enhanced risk for men to have prostate cancer [16]. Additionally, 
rs20576 variant carriers exhibited significantly enhanced colorec- 
tal cancer risk dependent on allele dose for female and also with 
advanced colorectal cancer stages [22]. However, some studies 
have also reported null association of this SNP with bladder [19], 
breast [10,42] and lung cancers [28,40]. In meta-analysis study, 
Chen et al also indicated increased risk of all types of cancer with 
AC and CC variants of rs20576 [51], but we failed to show any 
association with overall cancer risk and subgroups. 

We found the association of DR4 rs6557634 with GBC risk. 
Previously, Mittal RD et al., 2011 [19], also reported significant 
association of variant genotype of this SNP with bladder cancer 
risk. Though, Ulybina et al., 2009 [28] and Ulybina et al, 2011 
[42] failed to find any association with lung and breast cancer 
respectively. Our meta-analysis results showed a significant 
association of rs6557634 with overall cancer risk, GI cancers as 
well as in Caucasians. In-silico study using polyphen also showed 
this SNP to be possibly damaging. 

In our study, three haplotypes of DR4 C rs2 0575 

Ars20576A rs 6557634j G rs 20575A 1 - s 20576G rs 6557634 an d 

G rs 20575C rs 20576G rs 6557634 were associated with increased risk of 
GBC. The C rs 2057sA rs 20576A rs 6557634 haplotype persisted as risk 
factor even after gender or gallstone stratification of GBC cases. 
The C rs 20575-C rs 20576 diplotype was previously found to confer 2.4- 
fold risk of colorectal cancer by Frank et al., 2006 [22] and 3.5-fold 
risk of breast cancer by Frank et al., 2005 [10]. In another study, 
C rs 20575-G rs 6557634-C rs 20576 haplotype carriers had 1.8 folds 
increased risk in bladder cancer patients [19]. Thus, haplotypes 
seems to be better in predicting the association of DR4 with GBC 
risk rather than individual SNPs. 

FAS is a cell surface receptor found expressed in a variety of 
tissues and FASL is the natural ligand to FAS, whose expression is 
restricted to the activated T cells and natural killer cells [52] . The 
decreased expression of FAS may keep the transformed cells from 
elimination by antitumor immune response, called immune escape 
while the increased expression of FASL may increase the ability of 
tumor cells to kill FAS-sensitive lymphocytes, called immune 
counterattack [53]. Recently, up-regulation of FasL expression 
was reported in gallbladder carcinoma cells suggesting that FASL 
plays an important role in invasive depth, histological classification 
and metastasis of gallbladder carcinoma [54] . However, we failed 
to find any association of FAS (1377 G>A) and FASL (-844 C>T) 
polymorphisms with GBC susceptibility. Since these are receptor 
ligand system, an apoptotic cell death needs both normal FAS and 
normal FASL [14]. We have also looked for genotypic interactions 
between these two SNPs but we failed to find such an association 
(data not shown). Similarly, several other investigators have also 



failed to find an association of these SNP with cervical cancer [55], 
AML [56], lung cancer [57], gastric cancer [58] and NSCLC risk 
[27]. A meta-analysis by Zhang et al., 2009 [59] suggested that the 
FAS -1377G to A polymorphism may be a low-penetrance 
susceptibility marker of cancer. Another meta-analysis study by 
Zhang et al., 2009 [53] suggested a possible protective effect of 
FASL rs7631 10 T allele on cancer risk. These findings suggest that 
the promoter polymorphisms of the FAS and FASLG genes may 
not contribute to the etiology of cancer [58] and mutation of the 
primary structure of FAS or FASL might be one of the possible 
mechanisms that disrupt FAS-mediated apoptosis in tumor cells. 

FASL expression has been found to be induced by tobacco 
smoking [60]. DNA damage resulting from procarcinogenic 
compounds contained in tobacco smoke may increase DR4 
transcription and apoptotic cascade [61]. Considering the 
plausible combined effect of SNP and environmental exposure 
of various carcinogens like tobacco/ cigarette to modify individu- 
al's susceptibility to cancer, we have also looked for this in case 
only study. However, we did not found any association of the 
combined effect of tobacco and these SNPs in modulating the 
GBC susceptibility (data not shown). 

To the best of our knowledge, this is the first study investigate 
the effect of coding SNPs of DR4, FAS and FASL in predicting the 
GBC risk. The different effects of these aforesaid SNP on the 
various cancer risks may be partially ascribed to different 
molecular mechanisms underlying the pathogenesis of different 
cancers. Even in the same disease, different results may be due to 
difference in the genetic background in different ethnicities and/ or 
influence of different environmental factors. The allelic and 
genotype frequencies of these SNPs vary greatly in different ethnic 
groups. In addition, inadequate study design such as nonrandom 
sampling, limited sample size and the pitfalls arising from 
unknown confounders also need be considered. 

Conclusion 

DR4 haplotypes, especially C rs 20575A re 20576A rs6 557634 ! signifi- 
cantly increased the GBC risk and was unaffected by gender and 
gallstone status of patients. This haplotype may change the 
apoptotic signals, thus may modulate cancer susceptibility 
probably by promoting tumor cells survival and tumor growth 
rather than initiating tumor formation. It represents a late event 
and useful biomarker for GBC susceptibility. Though, it needs to 
re-confirmed in larger population based cohorts and validate in 
GBC patients of different ethnicities. In contrast, polymorphism of 
FAS and FASL system might not have any effect on influencing 
GBC susceptibility in North Indian population. Since a better 
understanding of the regulation of the signaling events and their 
perturbation in human cancers may lead to the identification of 
new molecular targets that can be exploited for therapeutic 
purposes. This strategy is expected to open new perspectives to 
target the death receptor pathway for cancer therapy. 

Supporting Information 

File SI Figure SI. Forest plot for rs20575. Figure S2. Forest 
plot for rs205756. Figure S3. Forest plot for rs6557634. Figure 
S4. Funnel plot for rs20575. Each dot represents an individual 
study for the indicated association. The horizontal lines represent 
CIs. Figure S5. Funnel plot for rs205756. Each dot represents an 
individual study for the indicated association. The horizontal lines 
represent CIs. Figure S6. Funnel plot for rs6557634. Each dot 
represents an individual study for the indicated association. The 
horizontal lines represent CIs. Table SI. Frequency distribution 
of DR4, FAS and FASL Gene polymorphism in GBC and HC 
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after subdividing on the basis of gender. Table S2. Frequency 
distribution of polymorphism in GBC and HC after subdividing 
on the basis of gallstone status. Table S3. Studies included in 
meta-analysis. Table S4: Functional information 
(ZIP) 
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